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Abstract 

We have calculated the two-loop strong interaction corrections to the neutralino 
pole masses in the DR scheme in the Minimal Supersymmetric Standard Model 
(MSSM). We have performed a detailed numerical analysis for a particular point 
in the parameter space and found corrections of a few tenths of a percent. We 
agree with previously derived analytic formulae for two-loop corrections to fermion 
masses. 



1 Introduction 



If supersymmetry (SUSY) as the most attractive extension of the Standard Model is real- 
ized at low energies, the next generation of high energy physics experiments at Tevatron, 
LHC and a future e + e~ linear collider will discover supersymmetric particles. Particu- 
larly at a linear collider, it will be possible to perform measurements with high precision 
[H [2] E] which allows us to test the underlying SUSY model. According to work done by 
the LHC/LC study group [3], the experimental accuracies for the determination of the 
neutralino masses of m^o and m^o are 50 MeV and 80 MeV, respectively, at the SPSla' 
benchmark point [3]. To match this accuracy it is indispensable to include higher order 
radiative corrections. 

In the framework of the Standard Model many calculations of quark pole masses were 
performed to two- loop [5]-[B] and three- loop order 0,01]] • Electroweak vector bosons were 
also studied in [H]-[IB] in this respect. Within the MSSM important results on quark 
self-energies were obtained in [H]-[2T]. Finally, in [22], [23] the gluino pole mass to two- loop 
order was calculated. 

General formulae for SUSY QCD corrections to fermion pole masses in the MSSM were 
derived in [22] • In this paper, we calculate the two-loop strong interaction corrections to 
the neutralinos within the MSSM in a completely independent way. We make a detailed 
numerical analysis for the benchmark point SPSla'. We agree with the analytic formulae 
given in [22] . 



2 Outline of the calculation 

The poles of the full renormalized propagator of a system of mixing Majorana fermions 
in the mass eigenbasis 

G m0j (k)- ^-. fc(1 + £, (fc2)) _ M t+s» 2 )|, 

in terms of the renormalized self-energies S^, ^m R , an d satisfy the complex pole 
mass equation 

= det( S - (M D - ± L m L ) ■ (1 + Elf)" 1 • (M^ - ■ (1 + i£ (2) 

Mr, = diag(m^o, m^o, m^o, m^o). The solution of this equation up to the considered order 
s P oie,i — TTT>ii — iTiTha = {M-dh + drrii) 2 determines the complex mass shift Srrii and thus 
the width Tj and the real part of the pole m^. In a perturbative expansion to 0(aots) for 
the case of neutralinos in the MSSM the shifts in terms of self-energies are written as 

brrti = 8m-/ + 5m^ 

= —(tg^i* + tTV + W L *M DU + ^M W ) 

—{t^i + t^ R V + t^mnu + t^ R :M DU ). (3) 
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where there is no sum on the index i. Note that there are no one- loop squared or one- 
loop differentiated contributions in eq. ([3]) at the two-loop level because as is absent at 
one-loop level. 
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Figure 2: Neutralino two-loop selfenergy diagrams with inner gluon line 
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Figure 3: Neutralino two-loop selfenergy diagrams with inner gluino line 

In Fig. Q] we show the one-loop diagrams. Note that diagrams with charge conjugated 
inner particles are not shown explicitly. We checked our analytic one-loop calculation 
against previous work p^l [25] in the on-shell scheme and found agreement. 

The nine different types of two-loop diagrams stem from SUSY QCD operator inser- 
tions into the quark/squark neutralino one- loop self-energy and can be grouped into those 
containing inner gluon (Fig. [2]), gluino (Fig. [3]) and squark lines (Fig. H]). The latter only 
contains the part of the 4-squark interaction proportional to as- Inserting (u, d) for every 
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Figure 4: Neutralino two-loop selfenergy diagrams with three inner squark lines 
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Figure 5: Neutralino two-loop self-energy diagrams with counter term insertions 

quark line in these diagrams gives a total number of 36 explicit two-loop diagrams to be 
calculated. Diagrams with one-loop counter term insertions (Fig. [3) involve 0{cts) mass 
counter terms for quarks and squarks as well as coupling constant counter terms stemming 
from the Yukawa part of the neutralino-quark-squark coupling and counter terms to the 
squark mixing matrix, see e.g. [19]. 

The diagrams and the amplitudes were generated in FeynArts 3.2 [26] and simplified 
analytically using FeynCalc 4.0.2 [27J. A small IR regulator mass parametrized artificial 
IR divergencies, and it was checked analytically that they cancel each other in the final 
result. According to the DR -scheme the UV- divergencies are regulated dimensionally 
by d = 4 — e, and the unphysical scalar mass parameter m t for the evanescent fields is 
absorbed according to [28] . The resulting integrals have a rather complex tensor structure 
which was reduced to a set of basis integrals using Tarcer [2S] . The resulting expressions 
are, however, too lengthy to be given here. The numerical analysis was then performed 
by implementing Tsil [3U] in a Fortran program which was auto-generated with the 
help of Mathematica. Apart from the check of gauge independence, we used the usual 
't Hooft Feynman R%=i gauge for the gluon field. 

3 Numerical results 

For the numerical analysis, we took the benchmark point SPSla' [I]. The SUSY pa- 
rameters at Q = 1 TeV are tan (3 = 10, M x = 103.209 GeV, M 2 = 193.295 GeV, 
M 3 = 572.328 GeV, // = 401.62 GeV, A t = -532.38 GeV, A b = -938.91 GeV, Mq 3 = 
470.91 GeV, = 385.32 GeV and = 501.37 GeV, for further details see again [3]. 

Fig. [6] shows the neutralino pole masses at SPSla' except for tan/3 which is varied 
between 3 and 30. At the SPSla' value tan/3 = 10 the corrections are larger than the 
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Figure 6: Neutralino pole masses m^o and rh^o as a funtion of tan (3 with experimen- 
tal accuracy at SPSla'. The red and the black line are the one- and two-loop results, 
respectively. 
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Figure 7: Neutralino pole masses as a function of \x for Xi ^nd \\ an< ^ the respective 
absolute corrections, all numbers in GeV. The red and the black line are the one- and 
two-loop results, respectively. Blue denotes the tree-level DR masses. 



experimental accuracy. 

The two-loop corrections to the masses of X3 an d x\ are of the order of 0.5 GeV lying 
far beyond the expected eperimental accuracy of 2 — 3 GeV [3] . 
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Figure 8: Left: relative corrections assuming gauge unification. Dashed and solid lines are 
one- and two-loop corrections, respectively. Right: Tree-level DR masses in this scenario. 
The black line is for Xi, blue for X2> re d f or X3 an d green for x\- Notice the level crossing 
between xfl an d x% at approx. 350 Ge V. 




Figure 9: Absolute neutralino mass shifts as a function of Mq, see text. All other pa- 
rameters from SPSla'. Left: absolute one-loop corrections. Right: absolute two-loop 
corrections. The black lines are for xi, blue for x^, re -d for x° and green for xfi- 



In Fig. [7] we vary fi from 300 GeV to 700 GeV. Since the tree-level neutralino masses 
strongly depend on this parameter we show the absolute one- and two-loop corrections in 
seperate plots. 

A rather different scenario is shown in Fig. [8] assuming gauge unification. The plot 
is over M2, \x is set to 200 GeV and all other values are taken from SPSla'. Taking 
over experimental accuracy from SPSla' which is 0.05% for m^o we again find important 
corrections in this scenario (one tick on the vertical axis is 0.2%). 

In Fig. [9] we show the one- and two-loop neutralino mass shifts as a function of the 
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Figure 10: Scale dependence of rh^o and m^o at SPSla'. The blue line is the tree-level 



DR input value, red and black are the pole masses at one- and two-loop level, respectively. 



third generation soft SUSY breaking masses Mq 
parameters are taken from SPSla'. 
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Mf) 3 . Again, all other 



4 Checks 

We have checked analytically the gauge-£ independence of the sum of those diagrams 
containing inner gluon lines in the 't Hooft gauge. The diagrams containing gluon lines 
provide e-scalar mass contributions which were treated separately and shown to drop out 
by redefining the soft bilinear squark breaking parameters [28] . 

The scaling behaviour given by the log[Q] and log 2 [Q] structure of our results was 
checked against the renormalization group prediction [21] analytically and numerically. 
For the latter, Fig. [TU] illustrates the improvement of the scale dependence for the pole 
masses of \\ and The two plots with only a red and a black line are zooms of the 
respective plots to their left in order to illustrate graphically the improvement when going 
from one- to two-loop level. For this plot we used the RGE subroutines of Spheno [32] 
which were extended to include also running standard model parameters. The remaining 
scale dependence of the neutralino pole masses indicate that further non-SUSY QCD 
corrections like Yukawa-interactions will give relevant contributions. 

The analytic RGE check is equivalent to the UV finiteness in the on-shell scheme. The 
trace of the square of the tree-level DR mass matrix, 



T DR = tr(Y^Y) DR = tr(MlM D ) DR = (M 2 + M 2 + 2(/i 2 + Aff )) 



DR 



(4) 



i=i 



is a sum of squared mass parameters and involves the whole set of two-loop self-energies 
when renormalized. In fact, it was shown that the sum of the squares of the pole masses 
is scale independent, 



o = A tr 

dt 



(M D (t) + 5m{X{t)))\M D {t) + 5m{X{t))) 



acts 



(5) 
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order by order in pertubation theory. X(t) abbreviates the set of DR input parameters 
X BR run to the scale Q. First, differentiating and expanding T(t) = tr(Y(t)'Y(t)) around 
t = log[Q/Qo] = with T(0) = T DR up to order aas and taas one has to evaluate 



dT(t) 
dt 




(6) 



using the respective one- and two-loop /3-functions /S^' 2 '* given in [3T] where sums over the 
indices of X and X' are understood implicitly. At the benchmark point SPSla' one has 
Qo = 1 TeV. The absence of strong interactions at the one-loop level is reflected by the 
fact that the as factor in taas on ly stems from P^l in eq. (jSJ). Second, the rhs of eq. (jfjj) 
has to be compensated by the Q-dependence of the mass shifts Srrii in eq. flSJ| evaluated 
with running parameters X(t). Note that in eq. (jSJ) the one- loop shift 5m^ has to be 
differentiated with respect to all ag-running parameters. Again, one-loop squared terms 
cannot contribute due to the absence of as at the one-loop level. 

As a last check, an independent C-program was written to evaluate the formulae given 
in [22] which agrees numerically with our auto-generated code up to the seventh relevant 
digit. 



5 Conclusions and Outlook 

Combining existing software packages [261 123 EH] we have calculated the two-loop 
corrections in 0(aas) to the neutralino pole masses in the MSSM in the DR scheme. We 
have shown that the results are important for matching theory and experiment at future 
colliders. Many checks have been succesful and the renormalization scale dependence 
is improved. However, it seems to be necessary to include further corrections, e.g. the 
leading two-loop Yukawa corrections 0(ah^,ahl). 
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